Stationary phase slip state in quasi-one-dimensional rings. 
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The nonuniform superconducting state in a ring in which the order parameter vanishing at one 
point is studied. This state is characterized by a jump of the phase by tt at the point where the 
order parameter becomes zero. In uniform rings such a state is a saddle-point state and consequently 
unstable. However, for non-uniform rings with e.g. variations of geometrical or physical parameters 
or with attached wires this state can be stabilized and may be realized experimentally. 
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In recent yearstl^H the existence of a single-connected 
state in a ring was proposed theoretically. Such a state 
implies that the relation between the phase cj) of the or- 
der parameter ?/; = /e"^, the current density j and the 
magnetic flux $ through the ring (which follows from the 
single valuedness of ■0 - see for example Ref.Q) 



(2b) 



—ds = 27rn - 
/ 



(1) 



is no longer valid (the flux is expressed in ^q/2tt, the 
current density in jo = c<i>o/47r^A^^, A is the London 
penetration length, ^ the coherence length and $0 is the 
quantum of magnetic flux). The reason is that the or- 
der parameter in such a single-connected state is zero at 
one point. Moreover, it was claimed that under certain 
conditions, i.e. radius of the ring less than ^ and the 
flux through the ring is about (n-l- l/2)$o this state inay 
become metastable in some range of magnetic fieldsH. 

In this paper we revisited this problem and we will 
show that a state where the order parameter vanishes in 
one point is still double-connected and that Eq. (1) is 
valid for such a state. We propose to call such a state a 
one-dimensional vortex state (ODV state), because like 
for an ordinary two-dimensional Abrikosov vortex, there 
is a jump in the phase of the order parameter of tt at 
the point where ^ = 0. In contrast to a two-dimensional 
Abrikosov vortex the ID vortex is an unstable structure 
in an uniform ring. However, if some inhomogeneities 
are present in the ring (defects, nonuniform thickness or 
nonuniform width of the ring, attached superconducting 
wires) the ODV structure can be stabilized and may con- 
sequently be realized experimentally. 

Consider a uniform ring with thickness d < A and 
width w < In addition, let the radius of the ring 
R be much larger than w. Under these conditions we can 
neglect the screening effects and the problem is reduced 
to a one-dimensional one. The distribution of the current 
density j and the order parameter -0 of the system at a 
temperature not far from the critical temperature Tc is 
described by the ID Ginzburg-Landau equation (plus the 
condition divj = 0) 



(Pf 



(2a) 



where the gauge-invariant momentum p = \7(j) — A is 
scaled in units of $o/(27r^), the length of the ring is 
L ~ I-kR and the circular coordinate s is in units of 
the coherence length ^. In these units, the magnetic field 
is scaled in units of the second critical field, iJc2, and the 
magnetic flux in <I>o/27r. 

The coupled Eqs. (2a) and (2b) have to be solved with 
the boundary conditipn ?/;(— L/2) ~ ^{Ljl). We use the 
method proposed inO (see alsoQ). These equations have 
the first integral 
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(3) 



From a formal point of view, Eq. (3) is nothing else then 
the condition of " energy" , {E) , conservation for some 
"particle" with "coordinate" / and "momentum" j. The 
role of " time" is playedQ by the circular coordinate s . In 
Fig. 1 we show the dependence of the "potential energy" 
of this system 
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for different values of j. Possible solutions of Eqs. (2a, b) 
are in the region where confined "trajectories" of our vir- 
tual particle exist. This is possible for currents < j < jc 
{jc — •\/4/27jo is the depairing current density ). Using 
Eq. (3) we immediately can write the solution of Eqs. 
(2a, b) for the ring 



V2s^ f 



dt 



-1 I t — to /<i — to 



\^t2 — to 



F sin 



^1 ^ ^0 ' V ^2 ~ ^0 



(5) 



where t{s) = f^{s), F{9,m) is the elliptic integral of the 
first kind and to < ti < t2 are the solutions of the cubic 
equation 



f' - 2f + AEt - 2f = 0. 



(6) 
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FIG. 1: Dependence of the "potential energy" U{f) for differ- 
ent values of the current j: 1- j=0.01, 2- j=0.05, 3- j=0.15, 
4- j^jc=^/I/27. 



Using the boundary condition for / we may conclude, 
that for a given current there exist a maximum of three 
sohitions. There are two uniform solutions (in the points 
of the minimum E = Emin and the maximum E — E^^ax 
of the "potential energy" (4)) and one nonuniform solu- 
tion with energy Emm < E < Emax which is defined by 
the equation 



L 



t2 — to 



-K 



ti - to 

^2 — to 



(7) 



where K{m) is the complete elliptic integral of the first 
kind. Numerical analysis of Eq. (7) for < j < jc shows 
that there is a minimal value for the ring length Lmm 
for which there exists a solution for this equation. When 
j jc, Ljnin ~^ oo and in the opposite limit j —> 
one can show that L„iin ^ tt. The latter corresponds 
to a ring radius R = 1/2 (or ^/2 in dimensional units). 
For a radius i? > 1/2 metastable states may exist and 
superconductivity is present for any value of tig magnetic 
flux (at least in our one-dimensional model)BlII. 

In principle, Eqs. (5-7) define the nonuniform solution 
of Eqs. (2a, b) for a ring. Unfortunately, even using the 
explicit Kordano expressions for the roots of Eq. (6) 
the results are rather complicated for arbitrary values of 
the current j. However a tractable analytical solution is 
possible if we consider the case of low currents j <C 1 for 
which the roots of Eq. (6) are simplified to 



to 



■ f/2E, 



1 - Vl - 4:E, 



t2-l + Vl - iE. 



(8c) 



After inserting these results into Eq. (7) we obtain the 
"energy". For L ^ tt the "energy" E ~ 1/4 is practically 



independent of L (e.g. for L = 8 the difference 1/4 — 
E{L ^ 8) is less than 0.002). Also in the hmit i - tt < 1 
we found that the " energy" E is independent of j 



(9) 



In the limit j ^ 1 we can also findEl the dependence of 
t(s) near the minimmn point of t{s) 



t{s) = 2f + L > TT, s < 1 



(10a) 



t{s)^f/2E + 2Es^, L-7r<l,s<l (10b) 

Using Eqs. (10a, b) it is easy to show, that in the limit 
j — > the gauge-invariant momentum is given by 



pis) = ilt{s) 



\J\ 



(vr + 27rm)(5(s), 



(11) 



where S{s) is the Dirac function and to is a integer. As 
a result the integral 



lim (p 4^ds 

3^0 J p 



lim 



(+^) 



: 



p{s)ds = 
±(7r + 27rTO), 



(12) 
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for arbitrary ring sizaa. Combining Eq. (12) with Eq. (1) 
it is easy to show that if $ = (fc + l/2)(f>o (fc is an integer) 
a solution of Eq. (2a) can be found which vanishes at one 

In previouE0fli studies of this state Eq. (11) was not 
taken into account (only the absolute value ofj-jbhe or- 
der parameter was found). Furthermore, in Ref.cl it was 
claimed that for small rings there is a magnetic field re- 
gion in which such a state exist and is metastable. But 
we find that in this region this state even does not exist 
(except in one point) because Eq. (1) is not fulfilled! 

The distribution of the order parameter as obtained 
in the present paper may also be applied to a ring in 
which some part of the ring consists of a normal metal. 
Then at the boundary between the normal metal and 
the superconductor the condition ^ ~ is fulfilled (if 
the normal part is longer than ^) and the distribution 
of the density |V'(s)P coincides with the one obtained in 
Ref.p (but the phase <j){s) will be different for our previ- 
ous ring geometry). BesidoSrwe cannot call such a state 
single-connected as in Ref.Ela'El because the relation (1) 
is valid even for this case. Therefore, it is better to call 
this state a one-dimensional vortex state because like the 
two-dimensional Abrikosov vortex there is a point where 
and the phase of the order parameter exhibits a 
jump equal to kn (as the Abrikosov vortex with orbital 
momentum 27rfc). A numerical calculation shows that, as 
in the case for an ordinary Abrikosov vortex in supercon- 
ductors where such a state is stable, the phase jump is 
always equal to tt, i.e. k = 1. 

Our numerical analysis of the time-dependent 
Ginzburg-Landau equations showed that this state is 
completely unstable for a uniform ring. However for rings 
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FIG. 2: Distribution of the absolute value (a) and phase (b) of 
the order parameter at different values of the magnetic flux: 
1- $=0.45$o, 2- $=0.49$o, 3- <&=0.498$o, 4- $=(0.5-0)$o, 
5- $=(0.5 + 0)$o, 6- $=0.502$o, 7- $=0.51$o, 8- $=0.55$o. 
Note that the order parameter for $ = <I>o/2 + a is the same 
as for "1> = <E>o/2 — a when a < $0/2. In the insets (a) fmin 
and (b) the phase difference Acp near the point s/L — 3/4 are 
shown as function of $. 



density will be equal to zero. It is interesting to note 
that the current density is also zero for $ — n^Q even in 
the case of a nonuniform ring, but for values of the flux 
vifhcre the ODV state does not exist. At these values of 
the magnetic flux the order parameter is uniform along 
the ring and the term 2Tm is completely compensated by 
the term $ in Eq. (1). 




1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
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FIG. 3: The current in the ring (with attached wire of length 
5^) as function of the flux. The different curves are for dif- 
ferent circumference of the ring: 1- L—1, 2- L—2, 3- L—3, 4- 
I/=3.5, 5- L=4. For I/=4 the ODV state does not exist in the 
ring and hysteresis appears. In the inset the dependence of 
dj/d^ on the flux $ is shown for rings with lengths L=l, 2, 
3, 3.5 (curves 1,2,3,4 respectively). Note that there are two 
maxima in the range (0, 1). With increasing length the two 
maxima merge into one and when one approaches the critical 
length dj/d^ diverges at $ = 'l>o/2. 



with nonuniform width(thickness) or attached wires a 
nonuniforra-diatribution of the order parameter becomes 
possibleEla'EHlEil and for small rings with L ^ n this state 
is realized in practice (besides for the case of a ring, a 
ODV-atate can also exist in the Wheatstone bridge - see 
RefsJlS-El). In Fig. 2 the distribution of the absolute 
value and the phase of the order parameter is shown for 
diffecent values of the magnetic field for a nonuniform 
ringlJ. We used the following model where the width of 
the ring was varied as 



w{s) 1 + lupsin 



2t:s 



(13) 



The parameters used are: wq ~ 0.1 and L — 3.25. It is 
seen that with increasing magnetic flux the magnitude of 
the order parameter / decreases in the thinner part of the 
ring (atjs/L — 3/4) and becomes zero when $/$o = 1/2. 
In Ref.B it was found that the ODV state is only possible 
at <& = {n + l/2)$o- The reason is now clear - only 
at this values of the magnetic flux the phase difference 
near the point where the order parameter is zero will 
be compensated by the term 27rn — <I> and the current 



Above we showed that the ODV state can be realized 
by varying the geometrical parameters of the ring. But 
there are two alternative approaches to realize the ODV 
state in the ring experimentally. First, it is possible to 
include an other phase in the ring. This leads to the ap- 
pearance of a weak link in the sample and if the radius 
of the ring is less than some critical value (about ^) it 
also leads to the existence of the ODV state in the ring 
at $ = (n -I- l/2)$o- We modelled this situation by in- 
troducing an additional term p{s)f in Eq. (2a), and, as 
an example, we took p(s) — —aS{s) and a = 1 (it should 
be noted that qualitatively the results do not depend on 
the specific value of a). We will not present the numer- 
ical results of the modified Ginzburg-Landau equations 
but we found that they are qualitatively similar to the 
behavior of / and (f> shown in Fig. 2. When approaching 
the flux <I> = (n -I- l/2)$o the order parameter reaches 
zero in the defect point and a jump in the phase equal 
to TT occurs. Secondly, such a ODV state should also ap- 
pear in | Sij ij ia| ll rings with attached wire(s). As was shown 
in Rcfs.lMjO an attached wire leads to a nonuniform dis- 
tribution of the order parameter in the ring. In Fig. 3 
the dependence of the current j in such a ring on the 
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FIG. 4: The Gibbs free energy and the magnetization (right 
lower inset) of the asymmetric ring (upper inset) as function 
of the applied magnetic field. 



flux through the ring is shown (see also Ref.0). As for 
nonuniform rings, where the parameters are chosen such 
that the OD\j-|State may exist, there is no hysteresis in 
such a systemEI. Small inhomogeneities in rings gives us 
the unique possibility to study states close to the ODV 
state because there is no sharp transition from the state 
with finite order parameter to a state with vanishing or- 
der parameter at one point (see Fig. 2). The latter is, in 
some sense, a 'frozen' slip phase state. A small variation 
of the flux through the ring near $ = $o/2 leads to a 
change of by 27r (from — vr to +7r with a concomitant 
change of the current from —0 to +0 - see Fig. 2(b)). 
As a result an additional phase circulation / \7(f>ds=2TT 
appears in the system, not because the magnetic flux <f> 
is changed, but because the term § j/f^ds changes from 
— TT to +TT (see Eq. (1)). When such a jump occurs V0 
does not change in the ring, except near the point where 
/ = 0, and as a result the current density in the system 
changes continuously. In the usual case transitions of the 
vorticity from a state with phase circulation 2Tm to a 
state with phase circulation 27r(n + 1) leads to a jump 
in the current everywhere in the ring because f "S/ (jyds 
changes by 27r and the order parameter and the current 
density are finite everywhere in the system. 

To supplement the above study we made also a nu- 
merical study of a nonuniform ring of flnite width and 
thickness by implementing our previous finite differ- 
ence solution of the coupled nonlinear Ginzburg-Landau 
equationalj. As an example we took the following pa- 
rameters: outer radius of the ring Ro = 2^, radius of 
the hole Ri = 1.5^, displacement of the hole from the 
center a = 0.4^, ring thickness d ~ 0.001^ and Ginzburg- 
Landau parameter k = 2. 

In Fig. 4 the dependence of the Gibbs free energy F 
and the magnetization M = —dF/dH of this system on 
the magnetic field are shown. As in the case of our one- 
dimensional ring these dependencies are reversible and 
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FIG. 5: The Gibbs free energy F{H) near the first maximum. 
In the inset the phase and the order parameter distribution 
at different values of the applied magnetic field (indicated by 
the squares on the F{H) curve) are shown. The phase was 
calculated along the outer perimeter of the ring. The open 
circle on F{H) indicates the position at which the vorticity 
increases from to 1. 



there are magnetic fields at which the magnetization is 
equal to zero (in the points of local maximum and mini- 
mum of the free energy). We found that the distribution 
of the order parameter and the phase in the ring (see 
inset of Fig. 5) is similar to analogical distributions for 
the above onc-dimcnsional ring (Fig. 2) at low magnetic 
fields. But different with the one-dimensional system the 
phase circulation increase of 27r (between the points (2) 
and (3) in Fig. 5) does not occur at the magnetic field 
value where the free energy has a local maximum (and 
zero magnetization). With decreasing width of the ring 
this point shifts towards the local maximum in the free 
energy. It is interesting to note that for the system corre- 
sponding to Fig. 4 the vortex enters through the thinnest 
part of the ring in the case of the first three maxima in 
the free energy while for the highest magnetic field max- 
ima (i.e. H/Hc.2 — 4.6) the vortex enters through the 
thickest part of the ring (see Fig. 6). Because the width 
of thicker part of the ring (= 0.9^) is considerable larger 
than the thinner part (= 0.1^) and is of order ^ the ODV 
for H/Hc2 — 4.6 becomes the usual two-dimensional vor- 
tex. As a result the order parameter is equal to zero only 
in one point along the radial coordinate and the circu- 
lation of the phase of the order parameter now is also a 
function of this coordinate (see Fig. 6). When the vortex 
enters/exits the ring at low magnetic fields there is also 
a slow dependence of the order parameter on the radial 
coordinate along the thinnest part of the ring. There- 
fore, we cair conclude that in a ring with finite width 
the one-dimensional vortex is transformed into the usual 
two-dimensional vortex. 

The process of the appearance of a stable vortex in a 
ring with finite width is similar to the creation of a phase 
slip in wires of finite width in the presence of a transport 
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FIG. 6: The same as Fig. 5 but now for the last maximum 
in Fig. 4. In the inset a contourplot of the order parameter 
is shown for different values of the magnetic field which are 
given by the squares in the main figure. L outer is the vorticity 
as calculated along the outer perimeter of the ring and Li„ner 
is the vorticity calculated along the perimeter of the hole. 



current. In the latter case the distribution of the gauge- 
invariant momentum is not exactly uniform over the wire 
width. The maximum value of p is obtained at the edges 
and as a result the order parameter vanishes first in these 
points. When a phase slip is created the distribution of 
the order parameter is not uniform over the width. With 
decreasing wire width this nonuniformity decreases but 
it will be uniform, strictly speaking, only in the limit 
w 0. With increasing wire width the phase slip trans- 
forms to the ordinary process of vortex/antivortex pair 
nucleating at the edges, penetrating deep into the su- 
perconductor and annihilating. In nonuniform rings the 
distribution of p over the width is nonuniform and in con- 
trast to wires with transport current it is nonsymmetric 
with respect to the ring width. Suppression of the order 
parameter first occurs only on the external side (when in- 
creasing the magnetic field) or on the internal side (when 
decreasing the magnetic field) . As in the case of an one- 
dimensional ring we may call this state a stable (" frozen" ) 
phase slip state if the width of the ring, where the vortex 
penetrates, is less than ^. 

The ODV state of a nonuniform ringjis very similar 
to the saddle-point state as found in Rof£3 in case of an 
uniform ring of finite width and in Ref.t£l for the case of 
a disk. Note that for an uniform ring with zero width 
the nonuniform solution Eq. (6) corresponds to a saddle- 
point point of F{tp). In Fig. 7 the gradual decrease of 
the hysteresis with increasing displacement a of the hole 
from the center of the ring is shown. It is seen that with 
increasing a the region where metastable states exist de- 
creases and ultimately vanishes for some critical displace- 
ment Ac. The saddle point states (only shown for a = 0, 
by the thin full curve) apppsaches the stable phase slip 
state when a = Qc- In Ref.tZI also nonuniform supercon- 
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C2 

FIG. 7: The Gibbs free energy F{H) at low magnetic fields 
for different displacements a of the hole from the center of 
the ring. The thin solid curve (for a — 0) corresponds to the 
saddle point state. 



ducting rings of finite width were investigated (but with 
larger ring radius). It was found that for low magnetic 
fields F{H) is irreversible but for high magnetic fields 
the dependence F{H) was reversible. In high fields the 
transition from a state with vorticity one to another state 
occured through the same scenario as discussed above. A 
similar state was also considered in Ref.Efl in the frame- 
work of the linearized Ginzburg-Landau equations for a 
nonuniform ring of finite width. They found that in such 
a system the vortex may be stable in some magnetic field 
region. In low magnetic field this vortex enter through 
the thinner part of the ring and at high magnetic fields 
through the thicker part. Our calculations generalizes 
this result to lower temperatures (i.e. where the non- 
linear Ginzburg-Landau equations have to be used) and 
with inclusion of the non-zero demagnetization factor of 
the ring. We found the dependence of the Gibbs free en- 
ergy (and magnetization) of this sample on the applied 
magnetic field. Our results also allowed to compare the 
results of the one-dimensional model with the full two- 
dimensional one. 

The ODV (or stable phase slip) state may be ob- 
served by magnetic experiments. Magnetic susceptibil- 
ity is proportionaB to dj/d^ and the ODV state exhibits 
some characteristic peculiarities as was shown in the inset 
of Fig. 3. Magnetization M is proportional to the cur- 
rent j and hence M{H) is reversible, for samples where 
the ODV state exists (see Fig. 3 and inset of Fig. 5). 
Furthermore M = at <i> ~ (n + l/2)<I)o- Alternatively, 
because &X ^ — {n + l/2)$o there is a point in the ring 
where the order parameter is equal to zero, this state 
may be^ound by transport measurements. For example 
in Ref.tj the dependence of the resistance of a hollow 
cylinder with radius of order ^ was studied at tempera- 
ture T < Tc far from Tc. At <i> = <i>o/2 the resistance 
(R) exhibited a maximum but the value of R was a fac- 
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tor three less than the normal state resistance i?„ . If the 
cylinder has a nonuniform thickness the ODV state can 
be realized as a stable state when $ = $o/2 and it will 
lead to a resistive (but superconducting) state even for 
very small currents and naturally the resistance of such 
a state will be less than i?„ as was found by Liu et altj. 

In conclusion, we studied the nonuniform state in a su- 
perconducting ring in which the order parameter vanishes 
at one point. It was shown that this state is characterized 
by a jump in the phase of the order parameter by tt near 
the point rp — 0. It allows, in correspondence with the 
ordinary two-dimensional Abrikosov vortex, to call such 
a state a one-dimensional vortex. This state is unstable 
in an uniform ring. In case of a nonuniform ring (with 
variations of the geometrical (d, w) or physical parame- 
ters (^, A) along the ring) or for a ring with an attached 
wire, this state may be realized in practise. For rings with 



non-zero width the ODV state transforms into the usual 
two-dimensional vortex. We also showed that this state 
is the remnant of the saddle point state connecting two 
superconducting states with different vorticity as found 
in a uniform ring. The latter state becomes stable in a 
nonuniform ring. 
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